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We consider triple Higgs boson production at a future 100 TeV proton-proton collider. We 
perform a survey of viable final states and compare and contrast triple production to Higgs boson 
pair production. Focussing on the hhh (bb) (bb) (jj) final state, we construct a baseline analysis 
for the Standard Model scenario and simple deformations, demonstrating that the process merits 
investigation in the high-luminosity phase of the future collider as a new probe of the self-coupling 
sector of the Higgs boson. 


I. MULTI-HIGGS BOSON PRODUCTION AT 
HADRON COLLIDERS 

With the exception of very few interactions, most 
of the terms that comprise the Standard Model (SM) 
Lagrangian have been measured or constrained, their 
strengths found to be suggestively close to the expected 
ones. An important category of interactions not directly 
observed are those of the the Higgs boson with itself. The 
so-called ‘self-couplings’ and their energy dependence are 
crucial in determining the stability of the vacuum. Cur¬ 
rent observations suggest that our Universe may be sit¬ 
ting at a metastable false vacuum HHsi and measure¬ 
ments of these couplings will illuminate this fact further. 

At colliders, these terms, i.e. those proportional to 
h being the Higgs boson scalar field, can be directly 
probed through the simultaneous production of (n — 1 ) 
Higgs bosons. Unfortunately, the production rates for 
processes with n > 3, i.e. more than one Higgs boson, are 
small, mainly due to the relatively large invariant mass of 
the final state system. In particular, at the Large Hadron 
Collider (LHC) with 14 TeV proton-proton centre-of- 
mass energy, gluon-fusion Higgs boson pair production is 
expected to have a cross section of ^40 fb isHn], whereas 
triple production is expected to have a rather dwarfish 
rate, with a cross section of 0(0.1 fb) [15]. Hence, even 
though there is optimism that Higgs boson pair produc¬ 
tion will provide important information and constraints 
through LHC measurements [TSHTS] , any direct measure¬ 
ment of SM-like triple Higgs boson production will be 
essentially impossible at the LHC, even at the end of 
the high-luminosity phase (HL-LHC) [76[[77]. However, 
with a significant increase in the collision energy, a Fu¬ 
ture Circular hadron-hadron Collider (FCC-hh), collid¬ 
ing protons at 100 TeV, stands a good chance at observ¬ 
ing and constraining the self-coupling of the Higgs bosons 
through Higgs boson pair production iiiniiisHHni, the 
cross section rising to ~ 1.6 pb | 8 l]. Additionally, at the 
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FCC-hh one may also get the chance to observe three on- 
shell Higgs bosons being produced, since the total cross 
section rises to ^5 fb [15]. The evaluation of this possi¬ 
bility is the main object of the present article. 

Concretely, the part of the Higgs boson potential which 
includes the self-interactions, may be written as: 

Vself = ^(l+C3)h3 + ^(l + d4)hy (1) 

where v ^ 246 GeV is the vacuum expectation value 
(vev), my 125 GeV is the measured Higgs boson mass 
and C 3 and 6/4 parametrize possible deviations from the 
standard model expectation (i.e. the SM is recovered for 
C3 = = 0). 

Figure shows some of the Feynman diagrams con¬ 
tributing to triple Higgs boson production. It is clear 
that the production cross section depends on both C 3 
and parameters. This should be contrasted to dou¬ 
ble Higgs boson production, which does not depend on 
(^ 4 . In Ref. [76] the dependence of the triple Higgs boson 
cross section on the parameters C 3 and ^4 was investi¬ 
gated at 14 TeV and 200 TeV proton-proton colliders 
for a Higgs boson mass my = 120 GeV. We produce an 
equivalent result for proton-proton collisions at 100 TeV, 
for my = 125 GeV, shown in Fig. The conclusions are 
similar to those drawn in US]: the cross section depen¬ 
dence on 6/4 is mild, the deviations due to ^4 = ±1 being 
at most ±20% for C 3 = 1. Hence modifications of the ^4 
coefficient itself will be very challenging to probe. This is 
also demonstrated in the contour plot of Fig.j^a), which 
shows the cross section normalised to the SM value, on 
the C 3 — ^4 parameter space. On this plane, one can ob¬ 
serve that the dependence along ^4 is much weaker than 
that along C 3 . 

In terms of constraining C3, triple Higgs boson pro¬ 
duction cannot be superior to double Higgs boson pro¬ 
duction due to its small production cross section. On 
the other hand, triple production would be the best pro¬ 
cess to constrain 6 / 4 , although, as we will demonstrate, 
even the FCC-hh with 30 ab“^ of integrated luminosity 
can only provide 0 ( 1 ) constraints on ^ 4 , because its de¬ 
pendency of the cross section is very modest. However, 
observing the triple Higgs boson production process is an 
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interesting task in its own right, and as will be seen, in¬ 
deed challenging at the FCC-hh. The goal of this article 
is to provide a first baseline study of Standard Model-like 
triple Higgs boson production via gluon fusion (ggF), at 
a future 100 TeV proton-proton collider. Furthermore, 
we investigate triple Higgs production in two scenarios 
where it is affected by new physics: (i) in the SM aug¬ 
mented by a single higher-dimensional operator in an ef¬ 
fective field theory approach and (ii) the generic case on 
the (c 3 — d 4 )-plane. 

The article is organised as follows: in Section EA] we 
investigate an explicit scenario that contains a single 
higher-dimensional operator. In Section [TT] we list, for fu¬ 
ture reference, the final states that could be interesting in 
the study of Higgs boson triple production. The Monte 
Carlo event generation, simulation of 6 -jet and photon 
tagging are described in Section [ml Differential distri¬ 
butions at parton level for triple Higgs boson production 
at 100 TeV, compared to those of Higgs boson pair pro¬ 
duction and the analysis of the channel ( 66 ) ( 66 ) ( 77 ) is 
described in Section |IV| We use this analysis to provide 
constraints in two scenarios. Finally, we provide discus¬ 
sion and conclusions in Section m 

A. The self-coupling in D = 6 EFT 

In the framework of the dimension -6 operator exten¬ 
sion to the Standard Model {D = 6 EFT), one can com¬ 
pare the sensitivity of multi-Higgs production to varia¬ 
tions of the operator Wilson coefficients [50]. Here we 
consider, as an illustrative example, a simplified mode 
with the assumption that the effect of all coefficients 
apart from a single one, originating from an operator of 
the form Oq ^ \H\^, where H is the Higgs doublet scalar 
before electroweak symmetry breaking: 

Veif = + MH\^ + Oe, Oe = (2) 

where jj? and A are the conventional parameters em¬ 
ployed in the SM potential for the Higgs doublet H. 

The changes in the quartic and the triple Higgs cou¬ 
plings, defined in Eq. are related via [50|:* 

C3 = C 6 , ^4 = 6 C 6 . (3) 

Due to the relation appearing in Eq.[^ the cross section 
for triple Higgs boson production is a quartic polynomial 
in C 6 , i.e. it contains terms up to Cg. Such terms come 
from squared matrix elements of diagrams containing two 
triple Higgs couplings, such as the one shown in Eig. Bd). 

In Fig. Ifb) we show the variation of the inclusive 
leading-order cross sections for ggE hh and hhh with re¬ 
spect to the SM (cg = 0). The fit as a function of cg for 


*Note that, in general, C3 and would be multiplied by j in 
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the two cases, at 100 TeV, is: 


^{pQ)hh 

= 0.22 X Cg 

cT{SM)hh 

0 


— 0.71 X Cg -|- 1.00, 

<^iC6)hhh 

= 0.03 X 4 

a{SM)hhh 

0 


T 0.03 X Cg -|- 0.43 X Cg 


- 1.31 X C6 + 1.00. 


The line = 6 C 3 is also shown as a dissection on the 
C 3 — plane in Eig. ia). 

II. TRIPLE HIGGS PRODUCTION FINAL 
STATES 

We list the dominant Higgs boson triple production fi¬ 
nal states, i.e. those that yield Vevents > 10 with 30ab“^ 
of integrated luminosity at a proton collider at 100 TeV 
centre-of-mass energy, in Table |l] 


hhh final state 

BR (%) 

a (ab) 

A^30ab-l 

{hb){hb){bh) 

19.21 

1110.338 

33310 

lbA{hh)lwWu) 

7.204 

416.41 

12492 

{bb){bb){Tf) 

6.312 

364.853 

10945 

{bb){Tf){WWu) 

1.578 

91.22 

2736 

{bb){bb){WW 2 e) 

0.976 

56.417 

1692 

{bb){WWu){WWu) 

0.901 

52.055 

1561 

(6b)(Tf)(Tr) 

0.691 

39.963 

1198 

{bb){bb)(ZZ 2 e) 

0.331 

19.131 

573 

{bb){WW2(){WWu) 

0.244 

14.105 

423 

(66)(66)(77) 

0.228 

13.162 

394 

{bb){Tf){WW2l) 

0.214 

12.359 

370 

{Tf){WWu){WWu) 

0.099 

5.702 

171 

{Tf){Tf){WWu) 

0.086 

4.996 

149 

{bb){ZZ2e)iWWu) 

0.083 

4.783 

143 

{bb){Tf){ZZ2i) 

0.073 

4.191 

125 

ibb)ijj){WWu) 

0.057 

3.291 

98 

(66)(Tf)(77) 

0.05 

2.883 

86 

{WWu){WWu){WWu) 

0.038 

2.169 

65 

{Tf){WW 2 e){WWie) 

0.027 

1.545 

46 

(T'f)(r'f)(r'f) 

0.025 

1.459 

43 

lbb)(WW2e)lwW2e) 

0.017 

0.956 

28 

{wW2t){WWu){WWu) 

0.015 

0.882 

26 

{bb){bb)(ZZu) 

0.012 

0.69 

20 

( tt ) ( tt ) (WW2£) 

0.012 

0.677 

20 

(bb)(ZZ2e)(WW2e) 

0.011 

0.648 

19 

(Tf)(ZZ2e)(WWu) 

0.009 

0.524 

15 

(bb)(jj)(WW2e) 

0.008 

0.446 

13 

(Tf)(jj)(WWu) 

0.006 

0.36 

10 


TABLE I: The list of channels with Aevents >10 with 30 ab“^ 
and their branching ratios (BR). The subscript denotes 
the number of leptons x in the final state, originating from 
the di-bosons. The cross section used for pp hh at 100 TeV 
is ctnlo = ctlo X 2.0 = 5.78 fb, where a A-factor K = 2.0 has 
been applied to obtain an estimate of the NLO cross section. 
The number of events has been rounded to the nearest integer. 










3 



FIG. 1: Example Feynman diagrams contributing to Higgs boson triple production via gluon fusion in the Standard Model. 
The vertices highlighted with a blobs indicate either triple (blue) or quartic (red) self-coupling contributions. 



FIG. 2: Total cross section ratio normalised to the Standard 
Model values for gluon-fusion-initiated triple Higgs produc¬ 
tion at 100 TeV obtained by varying the cs and parameters 
independently (see Eq. . The Higgs boson mass was fixed 
to iTih = 125 GeV. The SM cross section at leading order is 
~ 2.88 fb. The NNPDF23_iilo_as_0119 parton density function 
set was used. 


If we apply further requirements to the final states 
listed in Table [H 


• to possess greater than 100 events at 30 ab ^ of 
integrated luminosity, 


• and all gauge bosons fully decay to leptons, 

then we are left with the following interesting final states: 
{bb){bb){bb), {bb)(bb){Tf), {bb){bb){WW2i)^ {bb){Tf){rf), 
( 66 )( 66 )( 77 ), {bb){rr){WW 2 £). In particular, the ex¬ 
pected combined number of events in the multi- 6 -jet 
and multi-r final states is ^45000 over the lifetime of 
the FCC-hh, and will most likely provide valuable in¬ 
formation on the triple Higgs boson process. In the 
present study we focus on the rare but clean final state 
( 66 ) ( 66 ) ( 77 ). 


III. EVENT GENERATION AND DETECTOR 
SIMULATION 

A. Detector simulation 

In the hadron-level analysis that follows, we consider 
all particles within a pseudorapidity of |? 7 | < 5 and 
Pt > 400 MeV. We reconstruct jets using the anti-/ct 
algorithm available in the Fast Jet package [82l[83], with 
a radius parameter of R = 0.4. We only consider jets 
with Pt > 40 GeV within I 77 I < 3.0 in our analysis. We 
consider photons within \r]\ <3.5 and pt > 40 GeV and 
100% reconstruction efficiency. The jet-to-photon mis- 
identffication probability is taken to be = 10 “^, 

fiat over all momenta above the pr cut and over all pseu¬ 
dorapidities.'t' We also consider the mis-tagging of two 
light jets to bottom-quark-initiated jets with a flat prob¬ 
ability of 1 % for each mis-tag, corresponding to a flat 
6 -jet identification rate of 80% and demand that they lie 
within \r]\ < 3.0. We demand all photons to be isolated, 
an isolated photon having ^^PT,i less than 15% of its 
transverse momentum in a cone of AR = 0.2 around it. 
Finally, no detector-smearing effects have been consid¬ 
ered. 


B. Event generation 

Events for the hhh signal samples have been gen¬ 
erated via the loop-induced module of the MadGraph 
5/aMC@NL0 package [84H88] . The SM loop model present 
in MadGraph 5/aMC@NL0 was modified to allow for de¬ 
formations of the Higgs boson triple and quartic self¬ 
couplings away from the SM values. All tree-level and 
next-to-leading order (i.e. matched via the MC@NLO 
method [89]) background processes have been gener¬ 
ated using MadGraph 5/aMC@NL0, apart from the di- 
Higgs plus jets {hh + jets) background, which was simu¬ 
lated using HERWIG++ in conjunction with the OpenLoops 


iNote that the HL-LHC expectation has the approximate form 
= 0.0093 X e- 0 - 036 pTj/GeV For a PT 40 GeV, this 
gives approximately ~ 2 x 10“^. Thus, the value employed 

here is expected to be a reasonable approximation to future detec¬ 
tor performance. 
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FIG. 3: Total cross section ratios normalised to the Standard Model values for gluon-fusion-initiated multi-Higgs production 
at 100 TeV. The Higgs boson mass was fixed to rrih = 125 GeV. The SM cross section at leading order is ~ 2.88 fb. On the 
left-hand panel we show a contour plot of the variation of the cross section ratio with respect to the cs and parameters (see 
Eq.0). On the right-hand panel one can see the variation with respect to the SM in a theory where the SM is extended with 
a Oq ~ \H\^ operator as in Eq. for both Higgs boson pair production (hh) and Higgs boson triple production ihhh). For 
both calculations, the NNPDF23_nlo_as_0119 parton density function set was used. 


matrix-element generator isa EO]. The default par- 
ton density functions were used in each case: for the 
signal and tree-level backgrounds (including /i/i+jets) 
the NNPDF23_nlo_as_0119 set was used, whereas for the 
NLO samples the NNPDF23_nlo_as_0118_qed set was em¬ 
ployed [91]. 

Due to the large cross sections and high-multiplicity 
final states present at a 100 TeV collider, we only gener¬ 
ate the tree-level processes to include true photons and 
true 6 -quarks at parton level. This implies that light ex¬ 
tra jets for these processes will be generated by the par- 
ton shower, for which we employ the HERWIG++ general- 
purpose event generator [921495] .^ Inevitably this intro¬ 
duces an uncertainty to the results presented herein, ren¬ 
dering any observables related to these light jets leading- 
log accurate.^ We do not expect this, however, to alter 
the main conclusions of this first, baseline, study. Fur¬ 
thermore, generation-level cuts that anticipate the anal¬ 
ysis cuts at hadron level are imposed on the b quarks and 
the photons. In the case of decaying resonances (i.e. h 
and Z bosons) no cuts are imposed. The phase-space 
cuts applied on the samples 6666 , 66667 , 666677 , bbjj are 
shown in Table El 


^Simulation of hadronization and the underlying event were also 
included. Eg. No simulation of pile-up events was considered. 
§The /i/i+jets process is the only exception, with the first jet being 
leading-order accurate [32]. 


observable 

PS cut 


PT,b 

>35 GeV, at least one 

> 70 GeV 

M 

< 3.2 


PTn 

>35 GeV, at least one 

> 70 GeV 

\Vi\ 

< 3.5 


ARjj 

> 0.2 


TTL'y'y 

G [90,160] GeV 



TABLE H: The phase-space (PS) cuts imposed on the back¬ 
ground samples 6666 , 66667 , 666677 , 6677 . 


At this point one should stress that even though 
NLO event generation matched to the parton shower has 
been largely automated, NLO calculations for the high- 
multiplicity final states, particularly with many coloured 
particles and complicated phase space cuts, remain chal¬ 
lenging at present. We hence apply a conservatively large 
fiat iF-factor of K = 2.0 to all the processes calculated at 
tree level, as well as the hhh and /i/i+jets loop-induced 
processes. This is a crucial point that should be ad¬ 
dressed in future studies at higher-energy hadron collid¬ 
ers, as such final states will become increasingly common. 

The analysis of the signal and backgrounds gener¬ 
ated for the final state ( 66 ) ( 66 ) ( 77 ) is presented in sec¬ 
tion [TVBl 
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IV. ANALYSIS 


A. Differential distributions 

We investigate the shape of the differential distribu¬ 
tions in Higgs triple production in the Standard Model. 
Here we keep the Higgs bosons stable and include parton 
shower effects. We compare the shape of the hhh distri¬ 
butions to those coming from the more familiar case of 
Higgs boson pair production {hh) at 100 TeV. 

Figure Qa) shows the transverse momentum of any 
single Higgs boson either in hh or hhh production, 
Evidently, the transverse momentum of a Higgs boson 
in hhh is softer than that of hh, peaking at ^ 100 GeV 
instead of ^ 150 GeV. 

In Fig.Qb) we show the the spectrum of the transverse 
momentum of the Higgs boson “system”, i-e. the 

triplet of Higgs bosons in hhh, and the two Higgs bosons 
in hh. One can observe that the pT,h^ is harder in hhh 
than that of the pair in hh. 

We examine the distance between two Higgs bosons, 
Ai?(h, h), in hh and hhh production in Fig. |^c). In 
the case of triple production the distance is calculated 
between any two Higgs bosons. The Higgs bosons in hh 
are found to be more back-to-back than those in hhh, as 
expected. 

Finally, in Fig. |^d) we show the the invariant mass 
of all Higgs bosons in hh or hhh production, M^n. 
The invariant mass distribution in hhh peaks just above 
M /^3 ^ 600 GeV, whereas that in Higgs pair production, 
just above Mh 2 ^ 400 GeV. 

B. hhh^ (56)(66)(77) 

The hhh ( 66 ) ( 66 ) ( 77 ) process is expected to be rel¬ 
atively clean and simple to reconstruct.^ The excellent 
resolution of the di-photon invariant mass, that has con¬ 
tributed to the Higgs boson discovery at the LHG’s Run 
1 , can be exploited to facilitate background rejection. 

The present analysis follows a simple path, using the 
R = 0.4 anti-ht jets as described in Section |ml Note, 
however, that an analysis utilising the jet substructure 
of boosted Higgses to a bottom-anti-bottom pairs, e.g. 
as in [98], could assist in signal-background separation. 
We defer this task to future work. 

We ask for four 6 -jets, or light jets mis-identified as 
6 -jets, within I 77 I < 3.0, possessing transverse momenta 
Pt, 6 {i, 2 , 3 , 4 } ^ {80,50,40,40} GeV, where the subscripts 
1, 2, 3, 4 denote the first, second, third and fourth 
hardest 6 -jets respectively. We ask for two photons, 
or mis-identified jets as photons, within |? 7 | < 3.0 and 


^Note that this final state has been considered in EH , in the context 
of the two-Higgs doublet model hH —)► hhh final state. Here we 
consider the SM case. 


Pt, 7 {i 2 } > {70,40} GeV. Due to the fact that, for the 
majority of 6 -jets we cannot identify whether they orig¬ 
inated from a 6 -quark or an anti- 6 -quark, there exists a 
3-fold combinatorial ambiguity in combining 6 -jets into 
the two Higgs boson candidates. As a simple choice, we 
take the highest-pT ^-jet and pair it with the closest 6 - 
jet in AR = + A(/)^, and pair the other two re¬ 

maining 6 -jets together.il We thus construct the paired 
6 -jet invariant mass, respectively, and 

for which we demand G [100,160] GeV and 

^ciose ,2 ^ [90,170] GeV. The rather large mass win¬ 
dows are chosen to maintain high signal efficiency given 
the small initial cross section. Moreover, we construct 
the distance between the highest-pT ^-jet and the corre¬ 
sponding paired one, and impose ARf^^^'^ G [0.2,1.6].* ** 
For the photon pair, we simply construct the invariant 
mass and impose a strong window on the measured Higgs 
boson mass G [124,126] GeV.^^ We also restrict the 
distance between the two photons to AR^^ G [0.2,4.0]. 
We collect these selection cuts in Table iHll 


observable 

selection cut 

{1,2,3,4} 

> {80,50,40,40} GeV 

\Vb\ 

< 3.0 


e [100,160] GeV 

close,2 
^bb 

e [90,170] GeV 

A Dclose,! 

e [ 0 . 2 ,1.6] 

A oclose,2 

^^bb 

no cut 


> {70,40} GeV 

l»?7l 

< 3.5 

ARjj 

€ [ 0 . 2 ,4.0] 

TTi'y'y 

e [124,126] GeV 


TABLE III: The final selection cuts imposed in the analysis 
of the ( 66 ) ( 66 ) ( 77 ) final state. The observables are defined in 
the main text. 


We show a summary of the processes considered in the 
analysis in Table [IV| The most significant backgrounds 
in our set-up turn out to be the SM 666677 those 
coming from Higgs boson pair production in association 
with extra jets. Specifically, the latter emulates the signal 
well, as the di-photon mass window is expected to have 
similar efficiency to the signal. Moreover, as we have 
pointed out at the beginning of the section, the Higgs 
bosons in hh are harder on average than than those in 


II We have verified explicitly that an alternative method based on 
minimization of the squared sum of (mj,^ — m^) from each combi¬ 
nation yields results that differ by 0(1%) compared to the simpler 
AR method. 

**The distance between the other paired 6-jets was not found to have 
significant discriminatory power. 

’I”!'This cut implies that the di-photon resolution should be better 
than ~ 1 GeV at the FCC-hh. The current resolution at the LHC 
is 1-2 GeV, [99l [TOO] and thus it is not unreasonable to expect an 
improvement at the detectors of the future collider. 
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FIG. 4: Comparison of differential distributions for Higgs boson pair (hh) and triple production (hhh) in the Standard Model. 
Parton showering effects are included on top of leading-order matrix elements. Figure (a) shows the transverse momentum of 
any single Higgs boson, pT,h- In (b) we show the the spectrum of the transverse momentum of the Higgs boson “system”, 
i-o. the triplet of Higgs bosons in hhh, and the two Higgs bosons in hh. In (c) the distance between two Higgs bosons, 
AR{h, h), is examined and in (d) we show the the invariant mass of all Higgs bosons, Mh^. 


hhh, thus passing transverse momentum cuts easily. This 
background could be tackled in future studies via h ^ hh 
tagging using jet substructure techniques that exploit the 
decay versus the g ^ hh branching that produces the 
additional hb pair in /i/i+jets.^^ 


^^Note that the additional two 6-jets in /i/i+jets and hZ have been 
generated by gluon splitting into 66, performed by the shower 
Monte Carlo. 


C. Sensitivity in D = 6 EFT 

Despite the rather large backgrounds, a signal-to- 
background ratio of 0(1) can be obtained for the SM 
case. To summarise the results of the analysis, we present 
in the first two columns of Table [V| respectively, the num¬ 
ber of expected hhh events and the total expected num¬ 
ber of events, for the SM, as well as for the two simple 
deformations obtained by including the D = 6 operator 
Oq, with coefficient values cq = ±1. The third column 
of Table |V| indicates that, if one assumes that the SM 
is the underlying theory, then ce = ±1 can be excluded 
at 95% C.L. or better, using hhh (55)(55)( 77 ) at the 
‘high-luminosity’ phase of the FCC-hh. 

Furthermore, we show in Fig. the expected exclu¬ 
sion region on the cq coefficient, as well as the expected 
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process 

ctlo (fb) 

cnlo X BR X Vtag (ab) 

^analysis 

Arcuts 
^^30 ab-l 

hhh^ ( 66 ) ( 66 ) ( 77 ), SM 

2.89 

5.4 

0.06 

9.7 

hhh ( 66 ) ( 66 ) ( 77 ), cq = 1.0 

0.46 

0.9 

0.04 

1.1 

hhh ( 66 ) ( 66 ) ( 77 ), ce = — 1.0 

7.94 

15.0 

0.05 

22.5 

666677 

1.28 

1050 

2.6 X 10“^ 

8.2 

hZZ, (NLO) (ZZ ^ ( 66 )( 66 )) 

0.817 

0.8 

0.002 

<C 1 

hhZ, (NLO)(Z ^ ( 66 )) 

0.754 

0.8 

0.007 

< 1 

hZ, (NLO) {Z ( 66 )) 

8.019 X 10^ 

1129 

(9(10-®) 

< 1 

66667 + j^ts 

2.948 X 10^ 

2420 

( 940 -®) 

C(l) 

6666 + jets 

5.449 X 10^ 

4460 

( 940 -®) 

< 1 

6677 + jets 

98.7 

4.0 

(940"®) 

< 1 

hh + jets, SM 

275.0 

592.7 

7 X 10-'‘ 

12.4 

hh + jets , Ce = 1.0 

153.8 

331.5 

0.001 

9.9 

hh + jets , C 6 = — 1.0 

518.2 

1116.9 

4 X 10-* 

13.4 


TABLE IV: The processes considered in the analysis of the (56) ( 66 ) ( 77 ) final state. The parton-level cross section, including the 
cuts given in the main text is given (if any), the analysis efficiency and the expected number of events at 30 ab“^ are given. A 
flat iL-factor of iL = 2.0 has been applied to all tree-level processes (including /i/i+jets) as an estimate of the expected increase 
in cross section from LO to NLO. The hZZ^ hhZ and hZ processes have been produced at NLO and hence no X-factor is 
applied. Even though the hhZ process depends on ce, we only consider the SM case, as it was found to be negligible after cuts. 


number of events after cuts, at 30 ab“^. The theoretical 
uncertainty on the expected number of events for the hh 
and the /i/i+jets processes was taken to be 40% and un¬ 
correlated between the two. The analysis efficiencies for 
hhh and /i/i+jets were individually fitted using points in 
the region C 6 G [—3.0,4.0].^^ We assume that there is neg¬ 
ligible uncertainty on the ‘other’ backgrounds, which are 
taken to consist of the 666677 66667 +jets processes. 

By examining the central values of the the grey exclusion 
band, we can see that the regions cq < —0.7 and cq > 3.0, 
as well as the intermediate region cq G 1.0, ^ 1.7], are 
expected to be excluded at 95% C.L. {2a). Moreover, due 
to the fast-rising hhh cross section, as a function of the cq 
coefficient in this simple model, the 5cr-excluded region 
lies close to the 2a outer regions: cq < —1.4, cq > 3.5. 
Note that the analysis can be optimised for each value 
of Cq to obtain a higher significance, but in light of the 
many sources of uncertainties we do not pursue this here. 
Such optimisation could substantially alter the shape of 
the hhh and /i/i+jets curves in Fig. 


D. Sensitivity on the (c3 — d4)-plane 

Higgs boson triple production can be used to place con¬ 
straints on the (c 3 — (i 4 )-plane. This can subsequently be 
used to impose constraints on arbitrary relations between 
the triple and quartic coefficients in explicit models. We 
approximate the hhh signal efficiency over the whole 
plane by calculating its average value for C3 G [—3.0,4.0], 
<^4 = OC 3 , as obtained in the D = 6 EFT example. The 
analysis is used verbatim, without any modification of 


§§The fitting uncertainty is not shown in Fig. 



^6 

(a) 


FIG. 5: The expected exclusion significance on the cq coef¬ 
ficient (right vertical axis), assuming that the theoretical un¬ 
certainty on the expected number of hhh and hh+jets events 
is 40% for each process and uncorrelated between the two. 
The left vertical axis shows the expected number of events 
after cuts at 30 ab“^. The horizontal magenta dashed lines 
show the 2(7 and 5(7 exclusion points. 


cuts along the plane. The standard deviation on the ef¬ 
ficiency obtained this way was found to be ^ 20 % along 
this direction in the given interval. Considering the mag¬ 
nitude of the uncertainties on the signal and background 
predictions, we consider this to be adequate at present. 
For the /i/i+jets background we use the efficiency fit cal¬ 
culated for the D = 6 EFT case. We show the projected 
constraints on the (03 — (i 4 )-plane an integrated luminos¬ 
ity of 30 ab“^ in Fig. As a sanity check, we draw the 


















V. DISCUSSION AND CONCLUSIONS 



hhh 

total 

|iV(SM)-V(c6)| 

^JV(SM) 

SM 

9.7 

31.3 


C6 = 1.0 

1.1 

20.2 

- 2.0 

II 

1 

0 

22.5 

45.1 

2.5 


TABLE V: The number of events for an integrated luminosity 
of 30 ab“^ at 100 TeV, for the Standard Model and the the 
two simple deformations with Oq, with coefficient values ce = 
±1. The first and second columns show, respectively, the 
number of events for the hhh signal and the total expected 
number of events for all contributing processes: hhh, h/i+jets, 
bbbb'y'y (using 8.2 events) and bbbb'y-\-iets (using 1 event). The 
third column shows, approximately, the level (in number of 
standard deviations) at which the two hypotheses ce = ±1 can 
be excluded given that the standard model is the underlying 
theory. 


d 4 = 6 c 3 line and check that the outer 2 cr-region: cq ^ —2 
and cq > 3 approximately reproduces the D = 6 EFT 
result given the uncertainties. A few interesting observa¬ 
tions can be made. Firstly, the whole region C 3 < — 1 can 
be excluded at 5(7 irrespective of the value of <^4 using 
triple Higgs production. Moreover, if C 3 is constrained 
to lie near C 3 ^ 0 , then the weakest constraints on 
are obtained in all of the plane. On the other hand, if a 
non-zero value of C 3 is measured, e.g. C 3 ^ 4, then the 
constraint on can be quite stringent and in a region 
excluding d^ = 0, i.e. d^ G [ r\j r\j 8 ] at 5 ( 7 . 



FIG. 6 : The approximate expected 2a (blue) and 5(7 (red) ex¬ 
clusion regions on the c^ — d^ plane after 30 ab“^ of integrated 
luminosity, derived assuming a constant signal efficiency, cal¬ 
culated along the d^ — 603 line in C 3 G [—3.0,4.0]. 


Evidently, discovering Standard Model-like triple 
Higgs boson production will be a challenging task. Our 
analysis of the hhh (56) ( 66 ) ( 77 ) channel has demon¬ 
strated that the process merits serious investigation at a 
future collider running at 100 TeV proton-proton centre- 
of-mass energy. It is important at this point to emphasise 
the defining points and caveats that lead this phenomeno¬ 
logical analysis to this conclusion: 

• The detector of an FCC-hh needs to have excel¬ 
lent photon identification and resolution, so that a 
di-photon invariant mass window of width 2 GeV 
around the Higgs boson mass can imposed. As we 
already mentioned, the current resolution at the 
LHC is 1-2 GeV, [99l llQQj . Moreover, the pro¬ 
jections for photon identification efficiency at the 
high-luminosity LHC are at 0(80%) [IQlj . It is not 
unreasonable to expect an improvement in both of 
these parameters at the FCC-hh, to a resolution of 
< 1 GeV or photon identification of > 90%. 

• Tagging of 6 -jets should be extremely good, at least 
in the range of 70-80%, with excellent light jet re¬ 
jection of 0 ( 1 %) over a wide range of transverse 
momenta and pseudorapidities. Reducing the tag¬ 
ging probability from 80% to 70% would reduce the 
final number of events in ‘true’ 4-6-jet final states 
by about 40%. We note that the expected perfor¬ 
mance of the 6 -tagging algorithms for the LHC Run 
2 is already at this ballpark m- 

• Any analysis of triple Higgs production that in¬ 
cludes 66 pairs will also benefit from a very good 
forward coverage, allowing identification of 6 -jets 
up to pseudo-rapidities of I 77 I ^ 3.0. Good forward 
coverage for photons to I 77 I ^ 3.5 would also bene¬ 
fit the analysis. For example, the fraction of signal 
events with two 6 -jets falling in \r]i)\ G [2.5, 3.0] is 
^ 15% and the fraction of events with two photons 
falling in \r]j\ G [2.5, 3.5] is ^ 5%. These two are 
approximately uncorrelated, and thus an LHC-like 
coverage of < 2.5, \r]j\ < 2.5 would cause a 
^ 20 % reduction in signal efficiency compared to 
the analysis presented in this article. 

• Predictions of the triple Higgs boson production 
cross section, as for the case of double production, 
posses large theoretical uncertainties at present, 
due to the unknown higher-order corrections. The 
best available calculation includes only exact real 
emission diagrams in combination with ‘low-energy 
theorem’ results m- A full next-to-leading order 
calculation will reduce this and allow one to use the 
process to extract constraints on various models of 
new physics. 

• Crucially, the Monte Carlo event generation of mul¬ 
tiple coloured partons (4-6) at next-to-leading or- 
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der, with complicated phase-space cuts, matched 
to the part on shower, is essential. Technical im¬ 
provements in this direction, along with increase in 
computing power, will allow us to perform predic¬ 
tions with reduced theoretical uncertainties, as well 
as perform analyses of more hhh final states, such 
as those mentioned in Section (as well as other 
processes that involve multiple Higgs bosons). 

• Due to the aforementioned theoretical and techni¬ 
cal limitations, as well as the unknown characteris¬ 
tics of the future collider, we have not attempted to 
fully quantify the theoretical uncertainties perme¬ 
ating our results. We expect that future improve¬ 
ments in all of these aspects would allow one to 
obtain a more reliable quantitative result, includ¬ 
ing a reasonable expectation of uncertainty. 

We note here that our event selection is optimised for 
the assumed detector performance, and if some of these 
assumptions are changed, the event selection should also 
be changed to optimise the signal acceptance and back¬ 
ground rejection. Moreover in the scenario that the FCC- 
hh performance is substantially worse than what we have 
assumed, other channels could come into play, such as 
hhh {hh){hh){T^T~) or hhh (66)(t+t“)(t+t“). 

In conclusion, the study of triple Higgs production 
should be an important aspect of any future collider pro¬ 


gramme. It could provide complementary information on 
the nature of the Higgs boson and its role in electroweak 
symmetry breaking, as well as extensions of the Higgs bo¬ 
son sector beyond the standard model. This first baseline 
study resurrects this process and prompts further inves¬ 
tigation into how it can be put into use. 
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